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Abstract. We consider the local properties of the Yb 3+ ion in the crystal electric 
field in the Kondo lattice compounds YbRh 2 Si2 and YbIr 2 Si2- On this basis we have 
calculated the magnetic susceptibility taking into account the Kondo interaction in 
the simplest molecular field approximation. The resulting Curie- Weiss law and Van 
Vleck susceptibilities could be excellently fitted to experimental results in a wide 
temperature interval where thermodynamic and transport properties show non-Fcrmi- 
liquid behaviour for these materials. 
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1. Introduction 



Very peculiar magnetic, thermal, and transport properties of 4/ electron based heavy 
fermion systems are determined by the interplay of the strong repulsion of 4/ electrons 
on the rare-earth ion sites, their hybridization with wide band conduction electrons, and 
an influence of the crystalline electrical field. Consequences of the mentioned interplay 
for the electronic energy band structure near the Fermi-energy (Ep) were recently 
studied in YbRh 2 Si2 and YbIr 2 Si2 by angle-resolved photoemission and interpreted 
within the periodic Anderson model [TJ [2]. It was found that the hybridization of 
4/ electrons results in a rather flat 4/ band near Ep . Additionally, renormalization of 
the valence state leads to the formation of a heavy band that reveals strong 4f character 
close to Ep. Moreover, slow valence fluctuations of the Yb ion may occur between 
4/ 13 and closed 4/ 14 configurations with an averaged valence value of about +2.9 [3]. 
Evidently, these observations are consistent with a metallic behavior with very heavy 
charge carriers having properties of a Landau Fermi-liquid (LFL). At the same time the 
thermal, magnetic, and transport measurements show that the heavy fermions with a 
well defined Fermi-surface survive only at very low temperatures, coexisting with long- 
range antiferromagnetic (AF) order (Tjy = 70 mK in YbRh 2 Si2) which is suppressed at a 
magnetic quantum critical point (QCP) [4J by an external magnetic field (at H c = 600 G, 
Tjy — > K if c axis _L H ). With further increasing magnetic field at temperatures 
below a characteristic temperature T*, roughly proportional to H — H c , a crossover to 
LFL behavior is found. At temperatures above T/v and T*, but below the single ion 
Kondo temperature Tk, the properties of the discussed materials are quite unusual and 
display a non-Fermi liquid (NFL) behavior. The underlying fluctuations at the QCP are 
discussed to be locally critical, i.e. all the low-energy degrees of freedom have an atomic 
length scale jH [5] . One of the hallmarks of this local criticality is a generalized Curie- 
Weiss law like \ oc T~ a + const for the magnetic susceptibility with an exponent a < 1 
[6]. This type of behavior with a = 0.75 was found first in CeCus.gAuo.i [6]. The 
peculiarities of these properties are related to the competition between two interactions, 
originating both from the above mentioned hybridization: an exchange coupling of the 
local moments with the broadband conduction electrons (Kondo interaction) and an 
induced indirect RKKY interaction between the moments. The importance of the local 
properties in magnetic dynamics was mainly confirmed by the discovery of a strong and 
rather narrow electron paramagnetic resonance (EPR) in YbRh 2 Si2 and YbIr 2 Si2 below 
the Kondo temperature Tk = 25 K and Tk = 40 K, respectively (Tk revealed by specific 
heat data) [H El E]- This EPR signal was quite unexpected, since it was believed that 
the Yb 3+ magnetic moment at T < Tk should be screened by the conduction electrons 
and that the EPR linewidth should reach large values AH cx k B T K /giiB by approaching 
to Tk from above. Moreover, the main features of the observed EPR signal (anisotropy 
of the g factor and the EPR linewidth) reflect local properties of the Yb 3+ ion in the 
crystal electric field. The integrated intensity of the EPR line is proportional to the 
homogeneous static magnetic susceptibility. Having these experimentally confirmed 
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Table 1. Energies, wave functions and g factors of Yb 3+ ion in a tetragonal crystal 
field for T 7 and r| representations. 
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local properties at hand we assume entirely local properties of the Yb 3+ ion in the 
crystal electric field in order to theoretically investigate the energy spectra, g factors 
of the ground state, and the static magnetic susceptibility of YbRh 2 Si2 and YbIr 2 Si2 
compounds. 

2. Yb 3+ ion in tetragonal crystal field 

A free Yb 3+ ion has a 4/ 13 configuration with one term 2 F. As the spin-orbital coupling 
is much stronger than the crystal field in our compounds, the total momentum J is 
a good quantum number. The spin-orbital interaction splits the 2 F term into two 
multiplets: 2 F 7 / 2 with total momentum J = 7/2 and 2 F^j 2 with J = 5/2. Both are 
separated by about 1 eV [TU] and therefore we will consider only the ground multiplet 
2 Fj/2. The potential of the tetragonal crystal field for an ion can be written as 

V = aB° 2 0° 2 + (3{BlOl + S 4 4 4 4 ) + 7(5 6 °O 6 + 5 6 4 6 4 ). (1) 

To define energy levels and wave functions of the Yb 3+ ion we have to diagonalize 
the matrix of the operator ([T| on the states of the ground multiplet 2 F 7 / 2 - In 
(jl]) B\ are crystal field parameters, the operators 01(3) their matrix elements and 
a = 2/63, = -2/1155,7 = 4/27027 are given in ^Q\. 

As follows from the group theory, the two-valued irreducible representation D 7 / 2 of 
rotation group contains two two-dimensional irreducible representations Ttj and Tg of 
the double tetragonal group D 7 / 2 = 2Tj + 2Tq [10]. Therefore the states of Yb 3+ in a 
tetragonal field are four Kramers doublets and to diagonalize Hamiltonian Q we just 
need to diagonalize two two-dimensional matrices corresponding to the representations 
r^- and T 6 . Hence, the crystal field splits the lower 2 F 7 / 2 multiplet into four Kramers 
doublets with energies, wave functions and g factors as given in Table 1. In this table 
upper and lower signs correspond to left and right indexes; ], { correspond to the 
Kramers doublets effective spin projection up and down, c\ = cos(<^ c /2), c 2 = sin((/> c /2), 
a x = cos(y? a /2), a 2 = sin(y9 a /2), tany? a = A/A, tanyj c = C/C, -tt/2 < ip a , ^ c < n/2, 
and gj = 8/7 is the Lande g factor. We used parameters A, C, D, A, C which are defined 
by the crystal field parameters: 

A = 4B°/7 + 85^/77 + 80^/143, C = 4B°/2l + 405°/77 - 5605°/429, 
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D = 25°/21 - 645°/77 - 160^/429, 



(2) 



A = -8-^355^/385 + 80^^/3003, C = -8VSBt/77 - SOv^Sg/ 1287 - 

The Zeeman energy gjfisiiJ in the basis \ma) {m = 1...4, a =f , |) of each doublet could 
be represented by 



where H is the magnetic field, S is the effective spin operator with S = 1/2, hb is 
the Bohr magneton, gn and g± are g factors when the field is applied parallel and 
perpendicular to the c-axis, respectively (Table 1). 

As was mentioned above the main features of EPR signal observed in YbRh 2 Si2 and 
YbIr2Si2 reflect the local properties of Yb 3+ ion. The EPR signal in YbRli2Si2 and 
YbIr 2 Si2 is highly anisotropic [8j [IT] . The angular dependence of g factors in these 
compounds is well described by 



where 6 is the angle between magnetic field and crystal c-axis orientations, with 
|<7||| = 0.17, \g±\ = 3.56 for YbRh 2 Si 2 and \g»\ = 0.855, \g±\ = 3.36 for YbIr 2 Si 2 at 
T = 5 K. From neutron scattering experiments [12J [13] the intervals between the ground 
Kramers doublet and the exited energy levels amount to Ai = 17 meV, A2 = 25 meV, 
A 3 = 43 meV for YbRh 2 Si 2 and A 1 = 18 meV, A 2 = 25 meV, A 3 = 36 meV for 
YbIr2Si2- Unfortunately, these four independent values (three energy intervals and one 
parameter which define g» and g±, see Table 1) do not allow to determine five crystal 
field parameters unambiguously. 

Figure 1 represents the diagram of g factors together with the experimental points for 
the effective g factors of YbRh 2 Si2 and YbIr 2 Si2 (four points with different signs of 
g\\ and g±). The solid and dashed parts of the line gu + 2g± + 7gj = in Fig. 1 
correspond to the doublets E^Tq) and E^T^), and the solid and dashed parts of the 
ellipse (gn — gj) 2 /4 + gj_/3 = Ag 2 j correspond to the doublets £2^7) and EiiT^). It 
is evident from the proximity to the data (symbols in Fig. 1) that only the doublets 
£2^7) and E^Tq) could be considered for the ground state. The theoretical g values 
(i.e. the points on the ellipse and on the line nearest to the experimental points) and the 
corresponding values of ip c and ip a are given in Table 2. These results are qualitatively 
well consistent with the experimental ones being at the same time somewhat larger than 
the measured ones. Similar estimations for g factors were obtained recently in |T4j. 
A slight difference between experimental and our theoretical values can be explained 
mainly by taking into account the above mentioned the Kondo interaction, i.e. an 
exchange coupling between the 4/ electrons of the Yb 3+ ion and wide-band conduction 
electrons. This interaction becomes highly anisotropic after projection onto the ground 
Kramers doublet: 



^Zeeman — g\\^BH z S z + g±fiB(H x S x + H y S y ) 



(3) 




(4) 
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Figure 1. Diagram of g factors and experimental points at T = 5 K for effective g 
factors of YbRh2Si2 and YbIi'2Si2- 



Table 2. Optimal theoretical g factors of Yb 3+ ion in YbRh 2 Si2 and Ybrr 2 Si 2 . 



Compound 


Ground doublet, ip c , tp a 








511 


3-L 


YbRh 2 Si 2 


E 2 {T t 7 ),ip c = ±1.2660,- 
S 4 (r|), ip a = ±0.8206, - 


7T/2 < tp a 
-7T/2 < yj c 


< 
< 


it/2 
it/2 


-0.229 
-0.312 


±3.777 
-3.844 


YbIr 2 Si 2 


E 2 (r* 7 ), ip c = ±1.1003, - 
E 4 (r|), l p a = ±0.9952,- 


7T/2 < (/3 a 
-7T/2 < 99 c 


< 
< 


tt/2 
tt/2 


-0.929 
-0.940 


±3.529 
-3.530 



Here cr(rj) is the operator of the conduction electrons spin density, and J"' -1 are the 
exchange coupling integrals. This leads to the so-called Knight shift of the g factor. In 
the case of the anisotropic exchange interaction of the antiferromagnetic sign (j" ,J - < 0) 
this shift reduces the absolute value of the ionic g factor the same way as it happens in 
the isotropic case [15j [16] : 

9t = 9i±(l + V*0 > V = (6) 

y\\,±y<TiJ-B 

g?i ± is the ionic g factor, g a is the g factor of conduction electrons, Xa is the Pauli 
magnetic susceptibility, and X\\ t ± are molecular field constants. We also have performed 
an improvement of this simplest contribution of the Kondo interaction by the methods 
of a renormalization group analysis. Then, as will be published elsewhere, critical terms 
like In -1 (T/T^ 1 ) appear which reduce further the g factor values at low temperatures. 
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3. Static magnetic susceptibility 



The magnetization of the crystal is n(M), where n is the ion concen- 
tration, M = —gjHB>J is the ions magnetic moment operator, (M a ) = 
Tr(exp(— /37Y)M a )/Trexp(— /37^) is the mean value of the a component of magnetic 
moment, (3 = \/{k B T), k B is the Boltzmann constant, and T is the temperature. 
7i = V — MH is the Hamiltonian, where V is the crystal field potential ([!]) of the 
Yb 3+ ion, and the second term corresponds to the interaction of the magnetic moment 
with the magnetic field H. 
The magnetic susceptibility is defined as 



X«7 



n 



d(M a ) 



P 



n I (M a M y (X)) d\ } (7) 

H=0 Jo 

where (...) is calculated with V, M a {\) = exp(— \V)M a exp(AV). If we suppose the 
temperature to be such low that exp[— (3{E m — « for m ^ k we distinguish two 
different contributions in the susceptibility \ 
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Here a, a' =t,|- The first term, Xa^i corresponds to the Curie susceptibility 
proportional to inverse temperature, the second term, Xa-y > corresponds to the Van 
Vleck susceptibility which does not depend on temperature. \ka) indicates the states 
of the ground Kramers doublet. In the basis of the Kramers doublet states \ma): 
Xxz = Xyz = and Xxx = Xyy = XJ-> which is the evident result for tetragonal symmetry. 
We also introduce \zz = X\\> an d f° r the Curie constants C® z = CP, C® x = C® y = Cj_. 
As was shown above, the Kramers doublets ^(I^) and E^{T\) describe the ground state 
properties almost equally well. For these two cases the Curie and Van Vleck parts of the 
susceptibility could be expressed by the parameters Cj, a« and energy intervals between 
Kramers doublets. If the ground state is E2(T\) then 
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If the ground state is E^Tq) then 



xf = n/3(gjfjL B ) 2 [ 7y a 
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Table 3. Calculated Curie constant C° (10 6 m 3 mol X K) and Van Vleck susceptibility 
X vv (lO-^mor 1 ) for YbRh 2 Si 2 and YbIr 2 Si 2 . 



Compound 


Ground 
doublet 


cl 






4 V 


YbRh 2 Si 2 


Wr) 
Wa) 


16.8 
17.4 


0.087-0.202 
0.107-0.237 


0.062 
0.115 


0.09-0.227 
0.053-0.134 


YbIr 2 Si 2 


E2(n) 

Eitn) 


14.7 
14.7 


0.121-0.215 
0.127-0.224 


1.02 
1.04 


0.094-0.187 
0.083-0.166 



vv _ "i2n{gj\iB) 2 a\a\ 



X± V = n(gjfi B ) 2 



\a\a\ (^/7cia 2 — Vl5c 2 ai) 2 | (V7c 2 a2 + \/lhciai) 



+ 



E 3 — Ea 2(Ei — Ea) 2(E 2 — Ea) 



The calculated molar susceptibility values are given in Table 3. Here we used the 
g factors and parameters (p c and (p a from Table 2. In Table 3 the maximal and 
minimal possible values of X\\X. f° r different excited doublets sequences and uncertain 
parameters —tt/2 < (p a , (p c < tt/2 are given. We used for calculations experimental 
values of energy intervals Aj. It follows from our calculations that Curie- Weiss and Van 
Vleck susceptibilities play different roles in parallel and perpendicular orientations: for 
YbRh 2 Si2 the Curie constant in perpendicular orientation C° is at least two orders of 
magnitude larger than the Curie constant in parallel orientation Ci? whereas the Van 
Vleck susceptibility x± V nas ^ ne same order of magnitude as X» V ■ F° r YbIr 2 Si2 the 
situation is significantly different: C° ± and differ only by one order of magnitude 
whereas the Van Vleck part is almost the same as for YbRh 2 Si2. 

4. Comparison with the experimental data 

It is evident that the measured susceptibility includes both the Yb 3+ ions and the 
conduction electrons susceptibilities. In the molecular field approximation the Kondo 
interaction and RKKY interactions renormalize the total susceptibility [151 US E] : 

x,u + x ' = i-(^a+a l±Kl • *v = -r • < 12) 

where a\\,± are additional contributions to the molecular field from the RKKY interaction 

This renormalization leads to a Curie- Weiss law just as in the isotropic case [15] . The 
Pauli susceptibility is negligible. We neglected also the renormalization of the Van Vleck 
part of the susceptibility. Finally in the molecular field approximation we can write for 
the total magnetic susceptibility: 

* = TTtl + ^ (13) 
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C lu = Cl ± (l + 2\ lxX l) (14) 

and 6\\± independent on temperature. It is evident that we should expect C < C° if 
A < 0. 
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Figure 2. Static magnetic susceptibility of YbRh 2 Si2- Solid lines: fitting of 
susceptibility data with expression ( 13 1 and parameters shown in Table 4. Dotted 



and dashed lines: contributions of Curie- Weiss part of susceptibility in perpendicular 
(H = 11.2 kG) and in parallel orientations (H — 10 kG), correspondingly. Error bars 
indicate oricntational precision. Note the different scales for x± an d X\\- 



Figures 2 and 3 show the temperature dependence of susceptibility in YbRh 2 Si2 and 
YbIr 2 Si2 above 2K. It is easy to notice that the experimental data reflect the theoretically 
predicted tendencies according expression (13). Indeed, for YbRh 2 Si2 in perpendicular 
orientation the main role is played by the Curie- Weiss contribution to susceptibility, but 
in parallel orientation the susceptibility is almost temperature independent and the main 
contribution comes from the Van Vleck part. For YbIr 2 Si 2 the role of the Curie- Weiss 
susceptibility in parallel orientation is more important in comparison with YbRh 2 Si 2 . 
Table 4 presents the values of the fitting parameters C±, x± V '? @l and C\\, x\ V , 9\\ f° r 
YbRh 2 Si 2 and YbIr 2 Si 2 (as shown in Fig. 2 for YbRh 2 Si 2 Cn can be neglected within 
experimental error). As expected, the values of parameters C± and Cm are smaller than 
the calculated values C° and C». Indeed, as follows from (13), the renormalization of 
the susceptibility by the interaction with conduction electrons reduces the value of Curie 
constants because of the antiferromagnetic sign of the exchange integral (J" !± < 0) and, 
hence, A^j. < 0. From our fitting the Weiss temperatures are 9± = 5.43 K, 9u = 0.76 K 
for YbRh 2 Si 2 and 6± = 4.1 K, 0„ = 1.98 K for YbIr 2 Si 2 . 
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Figure 3. Static magnetic susceptibility of YbIi'2Si2 for H — 10 kG. Solid lines: fitting 



of susceptibility data with expression (13 1 and parameters shown in Table 4. Dotted 



and dashed lines: contributions of Curie- Weiss part of susceptibility in perpendicular 
and in parallel orientations, correspondingly. 



Table 4. Curie constants C (10 _6 m 3 mol _1 K), Van Vleck susceptibility x VV 
(10 -6 m 3 mol -1 ) and Weiss temperatures 9(K) for YbRli2Si2 and YbIr2Si2 from the 
data fits according Eq.((l3|) as shown in Figs. 2-4. 



Compound 


T(K) 


C x 


xl v 


0x 


G\\ 


4 V 


6\\ 


YbRh 2 Si 2 


2-30 
0.1-3.6 


10.89 
2.31 


0.064 
0.75 


5.43 
0.22 


0.02 


0.12 


0.76 


YbIr 2 Si 2 


2-30 


2.84 


0.18 


4.1 


1.04 


0.07 


1.98 



5. Discussion 

Our calculations of the magnetic susceptibility of YbRh 2 Si2 and YbIr 2 Si2 on the basis 
of an entirely local model of the Yb 3+ ion in the crystal electric field were stimulated 
by the following reasons. Firstly, the observed EPR signal reflects a number of features, 
which are very similar to those expected for Yb 3+ ions doped in non-conducting 
crystals (in particular, the local crystal field symmetry, the value of g factors, the 
temperature dependence of the EPR intensity, see [TO])- Even the EPR linewidth 
shows a temperature dependence that resembles the behaviour of Yb 3+ ions diluted 
in a conducting environment [15]. Secondly, an intensive experimental study of 
the NFL magnetic and thermal properties of these materials points out locally critical 
fluctuations. Thirdly, the calculations for the used local model could be performed in a 
straightforward and transparent way. 
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Figure 4. Fitting of the ac susceptibility data without additional dc field from [TH] by 
the power law p"5| with a — 0.6 (dotted line), a = 0.75 (dashed line) and by expression 
(p~3]) (solid line) with parameters shown in Table 4. 



Our major result is a remarkable agreement of our local approach for the static magnetic 
susceptibility with the temperature dependence of the experimental data. Therefore, in 
the considered region of temperatures (0.1 — 30 K) a ballistic motion of the 4/ electrons 
is practically absent, and they could be considered as quasi-localized. However, when 
approaching lower temperatures ferromagnetic quantum critical fluctuations dominate 
[U [18] and a locally quantum critical scenario may be applicable [5]. One of the 
hallmarks of this scenario is a generalized Curie- Weiss law which for a wavevector- 
dependent magnetic susceptibility can be written in the form 

C 

x(q,T) = TQ + g(q)Q (15) 

with an exponent a < 1 [5J In the case of YbRh 2 Si2 such a behaviour was revealed 
in the temperature region 0.3 < T < 10 K for q = 0, with a = 0.6 and 6 = 
[IS] . As shown in Fig. 4, we point out that a Curie- Weiss law together with a Van 
Vleck contribution convincingly describes the data for a wider temperature region in 



comparison with equation (15) and down to temperatures just above the AFM ordering 



temperature. However, for the temperature region shown in Fig. 4 the fitting parameters 
for the low temperature region are considerably changed, see Table 4. The reduction of 
the Curie constant C and Weiss temperature (8 = 0.22 K), as well as an increase of the 
temperature independent contribution can be related to the approach of the system to 
the LFL regime with a more ballistic motion of the 4/ electrons and indicate the Kondo 
effect in the magnetic susceptibility data. In this respect a Curie- Weiss description 
well within the Kondo regime, i.e. at T Tk, despite successful, may appear not 
appropriate. However, strong ferromagnetic correlations, as indicated, for instance, 
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by a large Sommerfeld- Wilson ratio for YbRli2Si2 [IE] and YbIr2Si2 dominate the 
magnetic susceptibility and may lead to this Curie- Weiss behavior. The reduction of 
the Curie constant can also be observed experimentally when comparing the magnetic 
susceptibility per Yb ion of YbRh 2 Si2 with Yi-^Yb^Pda (x = 0.6%) where the 4/ 
electrons are not hybridized with the conduction electrons [IS]. Interestingly, the Yb 3+ 
EPR intensity of the YPd3 Yb system compares well with the EPR intensity of YbRli2Si2 
P3J]. In respect of this, yet unexplained observation, it is worth to mention that in 
spite of the success of our entirely local approach for the static magnetic susceptibility 
of YbRh 2 Si2 this model is insufficient for a proper theoretical understanding of the 
dynamical susceptibilities as observed by EPR. We have found that this problem can 
be considered by taking into account a translational diffusion of 4/ electrons and 
their collective response together with wide-band conduction electrons to the resonant 
magnetic alternating field (the bottleneck regime). However, a discussion of this problem 
is beyond the scope of this paper and results will be published elsewhere. 
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